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Abstract 

Introduction: Expression of TXNDC5, wliicli is induced by liypoxia, stimulates cell proliferation and angiogenesis. 
Our previous study detected increased TXNDC5 expression in the synovial tissues of rheunnatoid arthritis (RA) 
patients using proteonnic nnethods. The current study investigated a pathogenic role for TXNDC5 in RA. 

Method: Expression of TXNDC5 in synovial membranes was quantitatively analyzed by immunohistochemistry, 
Western blotting and real-time polymerase chain reaction (PGR). Serum TXNDC5 levels and serum anti-TXNDC5 
antibody levels were determined using sandwich enzyme-linked immunosorbent assay (ELISA). A total of 96 single 
nucleotide polymorphisms (SNPs) in or near the TXNDC5 gene were genotyped using custom-designed lllumina 
96-SNP VeraCode microassay. Allele frequencies and genotype frequencies of SNPs were assessed using a case- 
control design in a cohort of 267 Chinese patients with RA, 51 patients with ankylosing spondylitis (AS) and 160 
healthy controls. Additional genotyping of 951 patients with RA and 898 healthy controls was performed for four 
SNPs (rs2277105, rs369086, rs443861 and rsl 1962800) using the TaqMan method. 

Results: Real-time PGR, Western blotting and immunohistochemistry detected significantly higher TXNDC5 
expression in the synovial tissues of RA patients compared to samples from patients with osteoarthritis (OA) or AS. 
ELISA detected significantly higher levels of TXNDC5 in the blood of RA patients compared to OA, AS and systemic 
lupus erythematosus patients, and healthy controls. ELISA did not detect significantly different levels of anti- 
TXNDC5 antibody in the blood of RA, OA and AS patients and healthy controls. A total of 9 SNPs (rs9505298, 
rs41 302895, rsl 225936, rsl 225938, rs372578, rs443861, rs408014, rs9392189 and rs2743992) showed significant 
association with RA, while 16 SNPs (rsl 0441 04, rsl 225937, rsl 225938, rs372578, rs89715, rs378963, rsl 225944, 
rsl 225947, rsl 238994, rs369086, rs408014, rs368074, rsl 225954, rsl 225955, rsl 3209404 and rs3812162) showed 
significant association with AS. Taqman SNP assay demonstrated that rs443861 has an association with RA, which 
correlates with the microassay results. 

Conclusions: TXNDC5 is up-regulated in synovial tissues of RA patients. TXNDC5 has a genetic effect on the risk of 
RA and AS. 



Introduction 

The thioredoxin domain, containing five (TXNDC5) pro- 
teins, also named ERp46, has a protein disulfide isomer- 
ase (PDI) domain that exhibits a high sequence similarity 
to thioredoxin, a catalyst of the rate limiting reaction of 
disulphide bond formation, isomerisation and reduction 
[1,2]. Yeast complementation tests showed that TXNDC5 
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can conduct PDI functions in vivo [3]. Indirect immuno- 
fluorescence microscopy and subcellular fractionation 
studies have shown that TXNDC5 is present both in the 
endoplasmic reticulum and the plasma membrane [4]. 
TXNDC5 is highly expressed in endothelial cells during 
hypoxic conditions, and plays important roles in anti- 
oxidative injury, anti-anoxia-induced apoptosis and the 
promotion of cell proliferation [1,2]. 

Abnormal proliferation of synovial fibroblasts and 
increased angiogenesis are pathological characteristics of 
rheumatoid arthritis (RA), an autoimmune disease that 
results in inflammation of the joints [5]. Using a proteo- 
mics approach, we detected increased TXNDC5 expression 
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in synovial tissues from RA patients [6] . Furthermore, we 
detected significantly elevated levels of TXNDC5 in the 
synovial fluid of patients with RA [6]. RA is thought to 
decrease the oxygen supply, leading to synovial hypoxia 
and hypoperfusion [7,8]. Hence, we believe that up-regula- 
tion of TXNDC5 may play an important role in the patho- 
genesis of RA in the hypoxic environment. 

In the current study, we quantitatively analyzed the 
expression of TXNDC5 in synovial tissues on both tran- 
scriptional and translational levels. We also examined 
TXNDC5 levels in the blood of RA patients using sand- 
wich ELISA. To determine genetic effects of TXNDC5 on 
RA, we conducted Illumina GoldenGate assays to identify 
potential associations between TXNDC5 polymorphisms 
and RA. SNPs, including tag SNPs, SNPs in promoter 
regions, SNPs in untranslational regions (UTRs), SNPs in 
exons and SNPs within proximity to exons of the 
TXNDC5 gene were genotyped in RA populations, and 
potential associations were determined by case-control 
study and haplotype analysis. 

Materials and methods 

Sample collection of synovial tissues and blood 

Synovial tissue samples were collected during knee joint 
replacement surgery from patients with RA {n = 10, 25 
female, 23 to 70 years old, mean 50) and patients with 
osteoarthritis (OA) {n = 10, 6 female, 41 to 77 years old, 
mean 60). Synovial tissue samples from patients with AS 
{n = 10, 3 female, 28 to 54 years old, mean 35) were col- 
lected during hip joint replacement surgery. The diagnosis 
of RA was made according to the criteria of the American 
College of Rheumatology. The patients with RA had dis- 
ease durations of 3-to-lO years and were classified as hav- 
ing erosive RA (Larsen class IV to V). They had high levels 
of C-reactive protein (30 to 100 mg/L, mean 24 mg/L), 
anti-CCP (300 to 3,000 U/ml) and RF (160 to 2,560 U/ml). 
AS patients had an average disease duration of seven years 
and were positive for HLA-B27 antigen. Their symptoms 
were consistent with the modified New York criteria for 
AS. Patients with AS and RA took disease-modifying anti- 
rheumatic drugs (DMARDs) before surgery. Patients with 
AS, RA and OA were also medicated with non-steroidal 
anti-inflammatory drugs (NSAIDs), which help reduce the 
pain and sweUing of the joints, and decrease stiffness. AU 
of AS and RA patients got treatment with DMARDs. 
Thus, the medical pretreatment does not influence the 
results and the experimental results are comparable. Addi- 
tional file 1 in the supplementary materials summarizes 
the epidemiological data. AU AS, RA and OA patients got 
treatment with NSAIDs. Synovial samples were dissected 
from connective tissues and immediately stored at -80°C 
until used. 

Peripheral blood samples were collected from patients 
with RA {n = 267, 183 female) and AS {n = 51, 10 female). 



RA patients had a mean age of 51.7 years, while AS 
patients had a mean age of 35.9 years. The diagnosis of 
RA and AS was conducted as described above. Patients 
were selected from the same population living in the Shan- 
dong area of Northern China. A total of 160 (58 female) 
healthy individuals with a mean age of 48.0 years were 
blood donors; they did not have any personal or family 
history of serious illness. Control individuals were fre- 
quency matched to the expected age distribution of the 
cases and were from the same geographical area. Blood 
samples were put into Monovette tubes containing 3.8% 
sodium citrate. 

Both patients and healthy controls gave their written 
consent to participate in the study and to allow their bio- 
logical samples to be genetically analyzed. The Ethical 
Committee of Shandong Academy of Medicinal Sciences 
approved this study. 

Western blot analysis 

Tissue samples weighing 200 (ig from RA, OA and AS 
patients were homogenized in Cell Lysis Solution 
(Sigma- Aldrich, St. Louis, MO, USA) and centrifuged at 
16,000 X g for five minutes at 4°C. Supernatants were col- 
lected after centrifugation, and protein concentrations 
were determined using the BCA Protein Assay Kit 
(Thermo Fisher Scientific, Rockford, IL, USA). Total pro- 
tein was separated by sodium dodecyl sulphate polyacry- 
lamide gel electrophoresis (SDS-PAGE) and trans-blotted 
onto nitrocellulose membranes (GE Healthcare, Piscat- 
away, NJ, USA). Western blot analysis was conducted 
using anti-TXNDC5 antibody (Abeam, Cambridge, 
Cambridgeshire, UK)) at a 2,000-fold dilution. The 
antibody was raised in goats using an oligopeptide 
(SLHRFVLSQAKDEL) against TXNDC5. All primary 
and secondary antibodies were diluted in 5% nonfat dry 
skim milk in TBST (Tris base 0.02 M, NaCl 0.137 M in 
distilled water (pH 7.6), containing 0.1% Tween-20). 
Immunoreactive signals were detected with alkaline 
phosphatase-conjugated secondary antibodies and visua- 
lized using a Western blotting luminol reagent (GE 
Healthcare). Western blot images were acquired on a 
Typhoon Trio (GE Healthcare). Quantification was con- 
ducted using ImageQuant 5.2 software. Another mem- 
brane prepared by the same protocol was probed with 
anti-GADPH antibody (Santa Cruz Biotechnology, Santa 
Cruz, CA, USA) to normalize sample loading. 

Immunohistochemistry 

Tissue sections of synovial tissues from RA, OA and AS 
patients were de-paraffinized and re-hydrated by standard 
procedures. Before the anti-TXNDC5 antibodies were 
applied, tissue sections were heated at 95°C for 10 minutes 
in citrate buffer solution (Sigma) for antigen recovery and 
then incubated with an endogenous peroxidase inhibitor 
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(Maixin-Bio, Fuzhou, Fujian, China) for 30 minutes at 
room temperature. After washing with PBS buffer (NaCl 
0.132 M, K2HP04 0.0066 M, KH2P04 0.001 5 M in dis- 
tilled water, pH 7.6), sections were incubated with antibo- 
dies directed against TXNDC5 (Abeam) overnight at 4°C. 
Immunoreactions were processed using the Ultrasensitive 
TM S-P Kit (Maixin-Bio) according to the manufacturer's 
instructions. Immunoreactive signals were visualized using 
DAB substrate, which stains the target protein yellow. Cell 
structures were counterstained with hematoxylin. 

In order to determine antibody specificity and opti- 
mize antibody dilution, the tissue samples were incu- 
bated (1) with goat pre-immune serum (Maixin-Bio, 
China) or (2) treated by the modification buffer without 
addition of antibody. 

Immunofluorescent labeling 

Tissue sections were processed as described above. After 
three washes with PBS buffer, tissue sections were treated 
with goat pre-immune serum (Maixin-Bio, China) for 
30 minutes to improve the specificity of the immunoreac- 
tion. Slides were incubated with anti-TXNDC5 antibody 
(Abeam) at 4°C for 12 h and then washed with PBS. 
TRITC 5-conjugated anti-goat IgG (Sigma-aldrich) was 
added to the slides, and slides were incubated for 40 min- 
utes at room temperature. Immunofluorescence was con- 
ducted with a Nikon 50i fluorescence microscope (Nikon, 
Shinjuku, Tokyo, Japan). To determine antibody specifici- 
ties and optimize antibody dilutions, a series of control 
slides were prepared as follows: primary antibodies only, 
secondary antibodies only and normal goat serum only. 
Expression levels of TXNDC5 were evaluated with Sim- 
plePCI (Hamamatsu Photonics, Sewickley, PA, USA), a 
semi-quantitative scoring system that analyzes the results 
of immunofluorescent labeling according to signal density. 

Real-time PGR 

Total RNA was isolated from the synovial tissues of RA, 
OA and AS patients using Trizol solution (Invitrogen Life 
Technologies, Carlsbad, California, USA) according to the 
manufacturer s protocol. Extracted total RNA was reverse- 
transcribed in a final volume of 10 \i\ using a RNA PCR 
Kit (TaKaRa, Katsushika, Tokyo, Japan). Real-time PCR 
reactions were conducted using the LightCycler 480 
Instrument (Roche Molecular Biochemicals, Basel, 
Switzerland) and performed according to the manufac- 
turer's protocol. Reactions were performed in a total 
volume of 10 ul, containing 1 ul of cDNA, 5 ul of SYBR 
Green Real-time PCR Master Mix (ToYoBo, Tokyo, 
Japan) and 1 ul of each primer. PCR amplification cycles 
were carried out as follows: 10 s at 95°C, 40 cycles of 5 s at 
95°C and 31 s at 60°C. For each sample, two reactions 
were performed at the same time. One reaction was per- 
formed to determine the mRNA level of the target gene. 



and the second was performed to determine level of P- 
actin. The experiment was performed in triplicate. PCR 
products were confirmed by melt curve analysis. Relative 
mRNA expression was calculated using the comparative 
threshold cycle (Ct) method according to the following 
formula: Ratio = 2-AACt = 2-ACt(sample), where ACt = 
Ct of target genes - Ct of endogenous control gene (P- 
actin). The relative target gene expression was normalized 
in comparison to P-actin mRNA levels. Primer sequences 
for the amplification of human TXNDC5 were as follows: 
forward primer for TXNDC5, 5'-GGGTCAAGATCGCCG 
AAGTA-3'; reverse primer for TXNDC5, 5'-GCCTCCA 
CTGTGCTCACTGA-3'; forward primer for human P- 
actin, 5'-TGGCACCCAGCACAATGAA-3'; and reverse 
primer for human P-actin, 5'-CTAAGTCATAGTCCGCC- 
TAGAAGCA-3'. Primer efficiency was determined by seri- 
ally diluting a standard RT reaction product. PCR 
efficiency was automatically calculated according to the 
dilution curve by the instrument software. Primer specifi- 
city was determined by both gel electrophoresis and melt 
curve analysis. Levels of TXNDC5 are expressed as the 
median and range. Statistical differences were assessed 
using the Mann- Whitney U-test; P < 0.05 was considered 
statistically significant. 

Sandwich ELISA detecting serum levels of TXNDC5 

Blood was collected from patients with RA {n = 96, 
75 females, 23 to 71 years old, mean 46), OA {n = 56, 
16 females, 50 to 86 years old, mean 62), AS {n = 56, 
19 females, 28 to 51 years old, mean 34) and systemic 
lupus erythematosus (SLE n = 56, 43 females, 23 to 
73 years old, mean 40) as well as healthy controls {n = 48, 
24 female, 20 to 40 years old, mean 31). Blood samples 
were collected using vacuum blood collection tubes. Fol- 
lowing centrifugation at 1,000 x g for 30 minutes, serum 
was collected and stored at -80°C until use. We raised 
antibodies in rabbits using an oligopeptide (RDGKKVD 
QYKGKRD) conjugated to keyhole limpet hemocyanin 
(KLH). The specificity of the antibody was confirmed by 
Western blot analysis using various recombinant proteins. 
The antibody was compared with the antibody made by 
Abeam, which showed similar results of immunohisto- 
chemistry and Western blotting. Rabbit antibody was 
diluted 5,000-fold in 0.05 M carbonate-bicarbonate buffer 
(pH 9.6) and used to coat 96-well ELISA microplates 
(Corning Life Science, Amsterdam, Netherlands) by over- 
night incubation at 4°C. After a brief wash with PBS con- 
taining 0.1% Tween-20 (PBST), plates were blocked with 
5% nonfat dry milk for one hour at room temperature. 
Next, blood samples were diluted 10-fold, and incubated 
in the plates for two hours at room temperature. After 
washing with PBST, goat anti-TXNDC5 antibody 
(Abeam), diluted 4,000-fold, was added to the plates and 
incubated for two hours at room temperature. Following a 
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washing step, a 15,000-fold dilution of anti goat IgG alka- 
line phosphatase-conjugated antibody (Sigma) was added, 
and plates were incubated for 30 minutes at room tem- 
perature. Following another PBST wash, plates were devel- 
oped by adding alkaline phosphatase yellow (pNPP) liquid 
substrate for ELISA (Sigma). Absorbance at 405 nm was 
measured using a plate reader (Synergy HT, Bio-Tek, 
Winooski, VT, USA). We repeated the ELISA three times 
and obtained the similar results. 

Sandwich ELISA has low inter-assay and intra-assay 
variability and provides more accurate results than 
direct ELISA in which patient sera were coated on the 
plate and were then detected using the antibody. 

ELISA detecting serum levels of anti TXNDC5 antibody 

Levels of anti TXNDC5 antibody were measured in the 
blood of patients with RA, OA, or AS {n = 50 for each dis- 
ease) as well as healthy controls {n = 50). One hundred 
microLs of SLHRFVLSQAKDEL (0.5 ug/ul), the oligopep- 
tide against TXNDC5, were coated onto 96-well ELISA 
microplates by overnight incubation at 4°C. After a brief 
wash with PBST, plates were blocked with 5% nonfat dry 
milk for one hour at room temperature. Serum samples, 
diluted 20-fold, were added and plates were incubated for 
two hours at 37°C. After washing with PBST, a 5,000-fold 
dilution of anti-human IgG alkaline phosphatase-conju- 
gated antibody (Sigma) was added, and plates were incu- 
bated for 30 minutes at room temperature. Following 
another PBST wash, plates were developed by adding the 
alkaline phosphatase yellow (pNPP) liquid substrate for 
ELISA (Sigma). Absorbance at 405 nm was measured 
using a plate reader. 

Genomic DNA extraction 

Genomic DNA was extracted from peripheral blood leu- 
kocytes using the DNA Blood Mini Kit from Qiagen (Hil- 
den, Germany) according to the manufacturer s guidelines. 
Briefly, 5 ml of blood was mixed with triton lysis buffer 
(0.32 M sucrose, 1% Triton X-100, 5 mM MgCl2, H2O, 
10 mM Tris-HCl, pH 7.5). Leukocytes were spun down 
and washed with H2O. Pellets were incubated with protei- 
nase K at 56°C and subsequently salted out at 4°C using a 
substrate NaCl solution. Precipitated proteins were 
removed by centrifugation. The DNA in the supernatants 
was precipitated with ethanol, and the resulting DNA pel- 
lets were dissolved in 400 [i\ H2O. 

SNPs selection 

lUumina GoldenGate assays were performed to genotype 
96 SNPs within or near the TXNDC5 gene in 267 RA 
patients, 51 AS patients and 160 healthy control indivi- 
duals from the Shandong area of North China. Tag SNPs, 
SNPs in untranslational region (UTR) and SNPs either in 
exons or in close proximity to exons of the gene encoding 



TXNDC5 were selected for genotyping. Tag SNPs were 
selected from HapMap data with a pair-wise r^ >0.8 and 
minor allele frequencies (MAF) over 0.05 [9,10]. Coding 
SNPs, SNPs near exons in 500 bp, SNPs in UTR and SNPs 
near the 5' and 3' ends of the gene were also selected. 
A total of 156 SNPs were candidates for lUumina's Gold- 
enGate design and were submitted to Illumina for a design 
score. The Illumina Assay Design Tool (Illumina, San 
Diego, CA, USA) filtered out SNPs not suitable for the 
Illumina platform, such as insertions/deletions, tri- and 
tetra-allelic SNPs, and SNPs that are not uniquely loca- 
lized. Finally, 96 SNPs with a design score of 1, spanning 
0.18 Mb of the chromosome were selected. These SNPs 
included 5 coding SNPs, 4 SNPs at the 3' UTR, 35 tag 
SNPs and 53 SNPs in introns or near the 5' end. The gene 
information of these SNPs is shown in Table 1. 



Table 1 Single nucleotide polymorphism (SNP) 
information 
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Table 1 Single nucleotide polymorphism (SNP) informa- 
tion (Continued) 
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Genotyping using microarray 

We performed genotyping using custom-designed Illu- 
mina 96-SNP VeraCode microarrays (lUumina). Genotyp- 
ing was completed by technique service in Dr. Zhang 
Feng's Laboratory of he Beijing Institute of Genomics. A 
BeadXpress Reader using Illumina VeraCode GoldenGate 
Assay Kit was used. A total of 500 ng of sample DNA was 
used per assay. Genotype clustering and calling were per- 
formed using BeadStudio software (Illumina). 

Genotyping using Taqman SNP assay 

Four tag SNPs, rs2277105, rs369086, rs443861 and 
rsl 1962800, were genotyped using TaqMan SNP genotyp- 
ing assays in a cohort of 950 patients with RA (693 female) 
and 900 healthy controls (630 female). RA patients had a 
mean age of 46.2 years and were from the Shandong area 
of Northern China. The diagnosis of RA was conducted as 
described above. Healthy individuals with a mean age of 
43.1 years were selected from the same geographical area. 
Assays were run on a LightCyclerH 480 Instrument 
(Roche) and evaluated according to the manufacturer's 
instructions. Reactions were carried out in a total volume 
of 10 [i\ using the following amplification protocol: dena- 
turation at 95°C for 10 minutes, followed by 40 cycles of 
denaturation at 92°C for 15 seconds and finishing with 
annealing and extension at 60°C for 1 minute. The geno- 
type of each sample was determined by measuring allele- 
specific fluorescence using SDS 2.3 software for allelic 
discrimination (Roche). Duplicate samples and negative 
controls were included to check the accuracy of 
genotyping. 

Statistical analysis 

Genotyping SNPs were analyzed for association by com- 
parison of the MAF in cases and controls. Associations of 
SNPs with RA and AS were evaluated using odds ratios 
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(OR) with 95% confidence intervals (CI). Fisher s exact 
test was used for comparison between categorical vari- 
ables. P-values less than 0.05 were considered statistically 
significant. After genotyping, SNP markers were evalu- 
ated for significant deviation from Hardy- Weinberg equi- 
librium. Calculations were performed using SHEsis and 
Haploview, two powerful web-based platforms for ana- 
lyses of linkage disequilibrium, haplotype construction 
and genetic association at polymorphism loci [11,12]. 

Results 

TXNDC5 expression in the synovial membranes of RA 
patients 

Immunohistochemistry analysis revealed significant 
TXCND5 expression in the thick lining layer and in 
many of the fibroblast-like cells of synovial membranes 
from RA patients {n = 10). Although detectable in the 
thin lining layer and some endothelial cells of small 
blood vessels, expression was very weak in the synovial 
membranes of OA patients {n = 10). In AS patients {n = 
10), TXNDC5 expression was relatively low in synovial 
membranes and was mainly limited in endothelial cells of 
small blood vessels. These observations were confirmed 
by immunofluorescent labeling. Results are shown in 
Figure lA. SimplePC (Hamamatsu Photonics, Sewickley, 
PA, USA), software designed to measure the signal den- 
sity of the expression in a semi-quantitative manner, 
detected significantly higher levels of TXNDC5 in syno- 
vial tissues from RA patients compared to OA and AS 
patients (Figure IB). 

Western blots revealed a protein with a molecular 
weight of 50 kDa. in each of the synovial tissues analyzed. 
Using GADPH as a reference, significantly increased 
TXNDC5 expression was detected in the synovial mem- 
branes of RA patients {n = 10), relative to the samples 
from OA {n = 10) and AS {n = 10) patients. These results 
were consistently observed in all of the synovial mem- 
branes examined (Figures 2 A, B). 

Transcription of TXNDC5 was quantified using real- 
time PCR. Similar to the Western blotting and immuno- 
labeling results, all RA samples {n = 10) exhibited a 
higher degree of TXNDC5 mRNA expression compared 
to the OA {n = 10) and AS {n = 10) samples (Figure 2C). 
TXNDC5 was expressed at a low level in all OA samples. 

TXNDC5 levels in blood samples from RA patients 

A sandwich ELISA was used to measure levels of 
TXNDC5 in the blood of RA patients with chronic 
inflammation. Levels of TXNDC5 were significantly 
increased in samples from RA patients compared to sam- 
ples from OA, AS and SLE patients. Serum TXNDC5 
expression in RA patients was also significantly elevated 
compared to healthy controls (Figure 3A). An ELISA was 
used to measure serum anti-TXNDC5 antibody levels of 



the patients. These were not significantly different from 
serum levels from RA, OA, and AS patients and healthy 
controls (Figure 3B). 

Genotyping of SNPs located in TXNDC5 

We genotyped 96 SNPs across the TXNDC5 gene from 
267 Han Chinese patients with RA, 51 patients and 160 
control individuals. All SNPs yielded genotype data, and 
the study sample success rate was 99.1%. Differences in 
allele frequencies and genotype frequencies between 
cases and controls were compared. Overall, nine SNPs 
(rs9505298, rs41302895, rsl225936, rsl225938, rs372578, 
rs443861, rs408014, rs9392189 and rs2743992) were 
found to be significantly associated with RA {P < 0.05). A 
total of 16 SNPs (rsl044104, rsl225937, rsl225938, 
rs372578, rs89715, rs378963, rsl225944, rsl225947, 
rsl238994, rs369086, rs408014, rs368074, rsl225954, 
rsl225955, rsl3209404 and rs3812162) were found to be 
significantly associated with AS {P < 0.05). Among the 
genotyped SNPs, three SNPs (rsl225938, rs372578 and 
rs408014) had significant association with both RA and 
AS. All SNPs retained in the analysis were in Hardy- 
Weinberg equilibrium {P > 0.05) in the overall samples. 
The allele and genotype frequencies of the associated 
SNPs between cases and controls are shown in Tables 2 
and 3. Other SNPs of the TXNDC5 gene polymorphisms 
did not disclose significant differences in allelic frequen- 
cies and genotype frequencies between the RA patients 
and controls or between AS patients and controls. 

Linkage disequilibrium (LD) analysis was performed 
within the tested SNPs. Pairwise D' values between all 
SNPs were calculated to determine the extent of LD. LD 
analysis defined eight blocks in TXNDC5 within the RA 
population. Rs372578, rs408014 and rs2743992, which 
showed strong association with RA, were in Blocks 3, 4 and 
8, respectively. LD analysis defined 10 blocks in TXNDC5 
within the studied AS population. Block 2 contained 
rs372578, rs89715, ra378963 and rsl225944, while Block 3 
contained rsl238994, rs369086, rs408014 rs368074, 
rsl225954, rsl225955, rsl3209404 and rs3812162, SNPs 
that showed strong association within AS patients. These 
results are shown in Figure 4A, B. 

In the RA population, haplotype analysis defined 27 
haplotypes (frequency > 1%) in the TXNDC5 gene by LD. 
Haplotype AA (frequency 79.3%) in Block 2, haplotype 
GAAG (frequency 56.4%) in Block 3, haplotypes 
GAGGGGA and AGCAAAC (frequencies 56.6% and 
23.1%, respectively) in Block 4 and haplotype AG (fre- 
quency 79.3%) in Block 8 provide significant evidence to be 
associated with RA risk {P = 0.0446, 0.0125, 0.0112, 0.0081 
and 0.0336, respectively). Haplotype analysis defined 40 
haplotypes (frequency > 1%) within the RA cohort by LD in 
the control population. Haplotypes AAAGAAG and 
GAAAGGA (frequencies 44.1% and 33.3%, respectively) in 
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Figure 1 Immunodetection of TXNDC5 in synovial membranes from patients with RA, OA and AS. (A) Immunolocalization of TXNDC5 in 
synovial membranes. The left lane indicates results of immunohistochemistry, and the right lane indicates results of immunofluorescent labeling. 
Original magnification: lOOx. Arrows indicate the upper layer of synovial membranes. (B) Semi-quantitative analysis of immunofluorescent signals 
of TXNDC5. TXNDC5 had significantly higher expression in the synovial tissue of RA patients compared to the synovial tissues of OA and AS 
patients. AS, ankylosing spondylitis; OA, osteoarthritis; RA, rheumatoid arthritis. 



Block 2, haplotypes AGGAGGGGA and CGAGCAAAC 
(frequencies 48.9% and 29.3%, respectively) in Block 3, hap- 
lotype GG (frequency 56.9) in Block 4, and haplotype AGA 
(frequency 4%) in Block 5 provided significant evidence to 



be associated with AS risk (P = 0.0198, 0.0043, 0.0044, 
0.0018, 0.0187 and 0.0053, respectively). The haplotype fre- 
quencies in a case-control cohort of patients with RA and 
AS are shown in Tables 4 and 5. The raw microarray data 
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Figure 2 Quantitative analysis of TXNDC5 expression. (A) TXNDC5 at molecular weight of 50 kDa was detected in synovial tissues of RA, OA 
and AS patients using Western blot analysis. Sample loading was normalized using GADPH at molecular weight of 37 kDa. (B) TXNDC5 
expression was semi-quantitatively analyzed by normalizing the signal density of TXNDC5 to that of GADPH. (C) TXNDC5 mRNA expression was 
measured in synovial tissues using real time PGR. The expression was normalized to that of p-actin. TXNDC5 had significantly higher expression 
in the synovial tissue of RA patients compared to the synovial tissues of OA and AS patients. AS, ankylosing spondylitis; OA, osteoarthritis; RA, 
rheumatoid arthritis. 
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Figure 3 Serum levels of TXNDC5 and anti-TXNDC5 antibody in 
patients with arthritic diseases and healthy controls. TXNDC5 
levels are represented by OD values of absorbance at 405 nm and are 
expressed as the mean ± standard error of the mean. (A) A sandwich 
ELISA detected increased level of TXNDC5 in blood samples from RA 
patients compared to samples from OA, AS and SLE patients, as well as 
from healthy controls. (B) An ELISA indicated that levels of anti- 
TXNDC5 antibodies were not significantly different among blood 
samples from RA, OA and AS patients and the healthy controls. AS, 
ankylosing spondylitis; OA, osteoarthritis; RA, rheumatoid arthritis. 



were available as an Additional file 2 to perform associa- 
tion, LD and haplotype analysis. 

We performed additional genotyping for four SNPs 
(rs2277105, rs369086, rs443861 and rsl 1962800) in an 
independent case-control study using the TaqMan 
method. The study was conducted within 951 patients 
with RA and 898 healthy controls. Allelic frequencies and 
gene frequencies of the four tag SNPs did not deviate 
from Hardy- Weinberg equilibrium in both case and the 
controls. Allelic frequency of the tag SNPs was compared 
between RA patients and controls. Among the polymorph- 
isms identified, the allele frequency and gene frequency for 
tag SNP rs443861 demonstrated statistically significant 
evidence for association with RA {P = 0.008320, and 
0.010110. This SNP was also determined to have signifi- 
cant association with RA by Illumina 96-SNP VeraCode 
microarray. The tag SNPs of rs2277105, rs369086 and 
rsl 1962800 did not disclose significant differences in alle- 
lic frequencies and gene frequencies between RA patients 
and controls (Table 6). 



Discussion 

In the present study, TXNDC5 expression was quantita- 
tively assessed both at the transcriptional level and trans- 
lational level. In comparison to synovial tissue samples 
from OA and AS patients, TXNDC5 expression was sig- 
nificantly increased in the synovial tissues of RA patients 
as determined by immunohistochemistry and Western 
blotting. Real time PGR also detected increased TXNDC5 
mRNA levels in the synovial membranes of RA patients. 
Furthermore, sandwich ELISA detected increased expres- 
sion of TXNDC5 in both the synovial fluid and blood of 
RA patients [6]. Taken together, these results confirm the 
increased expression of TXNDC5 in the synovium and 
blood of RA patients. In the present study, we did not 
detect increased levels of autoantibodies directed against 
TXNDC5 in the blood of RA patients, indicating that the 
over-expression of TXNDC5 does not directly cause an 
autoimmune response as an autoantigen like some citrul- 
linated proteins [13]. We processed Western blotting 
with protein extracted from the whole synovial tissue. 
The immunohistochemistry focuses on the expression of 
TXNDC5 in the lining area and the deep lining area of 
the synovial membrane. Immunofluorescent immunocy- 
tochemistry semi-quantified the expression level in one 
tissue region rather than the whole tissue. In addition, 
synovial tissues of RA and AS have significantly increased 
angiogenesis in which endothelial cells of blood vessels 
have strong expression of TXNDC5. Thus, it is possible 
the result of semi-quantification of immunofluorescent 
immunocytochemistry is a little different from the result 
of Western blotting. 

TXNDC5 expression is up-regulated by hypoxia and has 
a protective effect on endothelial cells by inducing folding 
and chaperone activity in hypoxia-induced anti-apoptotic 
molecules [1,2]. RA is thought to decrease the oxygen sup- 
ply, leading to synovial hypoxia and hypoperfusion [7,8]. 
Using co-immunoprecipitation followed by mass spectro- 
metry, Charlton et al found that TXNDC5 interacts with 
the N-terminal residues of AdipoRl. Further, transient 
knockdown of TXNDC5 in HeLa cells increased the levels 
of AdipoRl and AdipoR2, which correlated with the 
increased adiponectin-stimulated phosphorylation of 
AMPK. However, adiponectin-stimulated phosphorylation 
of p38MAPK was reduced following TXNDC5 knockdown 
[4]. Recent reports indicate that AdipoRl and AdipoR2 
mediate the insulin-sensitizing adipokine adiponectin. RA 
is associated with the increased production of adipokines, 
cytokine-like mediators that are produced mainly in adi- 
pose tissue and synovial cells [14]. Frommer et al demon- 
strated that adiponectin was present in inflamed synovium 
at sites of cartilage invasion in lymphocyte infiltrates and 
in perivascular areas. Adiponectin stimulates synovial 
fibroblasts to secrete chemokines, proinflammatory cyto- 
kines, prostaglandin synthases, growth factors and factors 
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Table 2 Allele and genotype frequencies in a case-control cohort of patients with RA 
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160 controls and 266 cases were observed; ""represents tag SNP. 

95% CI, 95% confidence Interval; RA, rheumatoid arthritis; SNP, single nucleotide polymorphism. 



for bone metabolism and matrix remodelling. This adipo- 
nectin-mediated effect was p38 MAPK and protein kinase 
C dependent. Adiponectin promotes inflammation 
through cytokine and chemokine production that attracts 



inflammatory and pro-destructive cells to the synovium, 
which, in turn promotes matrix destruction at sites of car- 
tilage invasion [15]. Choi et al. reported that adiponectin 
might contribute to synovitis and joint destruction in RA 
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Table 3 Allele and genotype frequencies in a case-control cohort of patients with AS 

dbSNP identity Allele/Genotype Numbers of patients Numbers of controls (%) Fisher's P-value Odds Ratio (95% CI) 



with RA (%) 


rs 1044 104 


A 


49 (0.480) 


103 (0.322) 


0.0037 


1 .947793 (.237327 to 3.066204) 




G 


53 (0.520) 


217 (0.678) 




1.947793 (.237327 to 3.066204) 




AA 


13 (0.255) 


20 (0.125) 


0.022555 






AG 


23 (0.451) 


63 (0.394) 








GG 


15 (0.294) 


77 (0.481) 






rsl 225937 


A 


67 (0.657) 


252 (0.787) 


0.007503 


0.516553 (0.316863 to 0.842090) 




G 


35 (0.343) 


68 (0.212) 




0.516553 (0.316863 to 0.842090) 




AA 


22 (0.431) 


99 (0.619) 


0.027945 






AG 


23 (0.451) 


54 (0.338) 








GG 


6 (0.118) 


7 (0.044) 






rsl 225938 


A 


46 (0.451) 


190 (0.594) 


0.011468 


0.562030 (0.358612 to 0.880834) 




G 


56 (0.549) 


130 (0.406) 




0.562030 (0.358612 to 0.880834) 




AA 


11 (0.216) 


52 (0.325) 


0.013973 






AG 


24 (0.471) 


86 (0.537) 








GG 


16 (0.314) 


22 (0.138) 






rs372578 


A 


50 (0.490) 


109 (0.341) 


0.006659 


1.861327 (1.184652 to 2.924518) 




G 


52 (0.510) 


211 (0.659) 




1.861327 (1.184652 to 2.924518) 




AA 


14 (0.275) 


21 (0.131) 


0.029107 






AG 


22 (0.431) 


67 (0.419) 








GG 


15 (0.294) 


72 (0.450) 






rs89715 


A 


52 (0.510) 


212 (0.662) 


0.005544 


0.529811 (0.337112 to 0.832662) 




G 


50 (0.490) 


108 (0.338) 




0.529811 (0.337112 to 0.832662) 




AA 


15 (0.294) 


73 (0.456) 


0.026647 






AG 


22 (0.431) 


66 (0.412) 








GG 


14 (0.275) 


21 (0.131) 






rs378963 


A 


63 (0.643) 


253 (0.791) 


0.002897 


0.476680 (0.291090 to 0.780596) 




G 


35 (0.357) 


67 (0.209) 




0.476680 (0.291090 to 0.780596) 




AA 


20 (0.408) 


100 (0.625) 


0.012477 






AG 


23 (0.469) 


53 (0.331) 








GG 


6 (0.122) 


7 (0.044) 






rsl 225944 


A 


35 (0.357) 


66 (0.206) 


0.002276 


2.138047 (1.304553 to 3.504071) 




G 


63 (0.643) 


254 (0.794) 




2.138047 (1.304553 to 3.504071) 




AA 


6 (0.122) 


7 (0.044) 


0.010406 






AG 


23 (0.469) 


52 (0.325) 








GG 


20 (0.408) 


101 (0.631) 






rsl 225947 


A 


51 (0.500) 


108 (0.338) 


0.0032 


1.962963 (1.249087 to 3.084832) 




C 


51 (0.500) 


212 (0.662) 




1.962963 (1.249087 to 3.084832) 




AA 


13 (0.255) 


21 (0.131) 


0.017854 






AC 


25 (0.490) 


66 (0.412) 








CC 


13 (0.255) 


73 (0.456) 






rsl 238994 


A 


52 (0.510) 


210 (0.656) 


0.007966 


0.544762 (0.346807 to 0.855709) 




C 


50 (0.490) 


110 (0.344) 




0.544762 (0.346807 to 0.855709) 




AA 


15 (0.294) 


71 (0.444) 


0.031615 






AC 


22 (0.431) 


68 (0.425) 








CC 


14 (0.275) 


21 (0.131) 






rs369086* 


A 


35 (0.343) 


66 (0.208) 


0.005322 


1.994573 (1.221370 to 3.257261) 




G 


67 (0.657) 


252 (0.792) 




1.994573 (1.221370 to 3.257261) 




AA 


6 (0.118) 


7 (0.044) 


0.022037 






AG 


23 (0.451) 


52 (0.327) 








GG 


22 (0.431) 


100 (0.629) 






rs408014 


A 


50 (0.500) 


211 (0.659) 


0.004146 


0.516588 (0.327716 to 0.814312) 
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Table 3 Allele and genotype frequencies in a case-control cohort of patients with AS (Continued) 
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70 (0.700) 


261 (0.816) 




1.895884 (1.135508 to 3.165434) 




AA 


4 (0.080) 


7 (0.044) 


0.043599 






AG 


22 (0.440) 


45 (0.281) 








GG 


24 (0.480) 


108 (0.675) 






rs3812162^ 


A 


70 (0.700) 


268 (0.838) 


0.002475 


0.452736 (0.268969 to 0.762060) 




C 


30 (0.300) 


52 (0.163) 




0.452736 (0.268969 to 0.762060) 




AA 


23 (0.460) 


114 (0.713) 


0.004576 






AC 


24 (0.480) 


40 (0.250) 








CC 


3 (0.060) 


6 (0.037) 







160 controls & 266 cases were observed; ^represents tag SNP. 95% CI, 95% confidence Interval; AS, ankylosing spondylitis; SNP, single nucleotide polymorphism. 



by stimulating vascular endothelial growth factor, matrix 
metalloproteinase-1, and matrix metalloproteinase-13 
expression in fibroblast-like synoviocytes [16]. Addition- 
ally, Tian et al also reported that increased PDI activity in 
myocardial endothelial cells in mice stimulates angiogen- 
esis under hypoxia condition [17]. These results support 
the possibility that the increase of TXNDC5 expression in 
the synovial tissues of RA patients stimulates the synovial 
ocular pannus, pro-inflammation and bone degradation. 
However, the detailed mechanism requires further 
investigation. 

TXNDC5 is a newly identified member of this protein 
family. TXNDC5 has been genetically mapped to chro- 
mosome 6p24.3. The gene encoding TXNDC5 is 
approximately 845.2 k bp, and it is divided into 13 
exons. The present study genotyped 96 SNPs flanking 
the TXNDC5 gene through Illumina GoldenGate assays. 
Further, the study also genotyped four tag SNPs in the 
TXNDC5 gene using the Taqman method to confirm 
association to RA in a large number of samples. Both 
methods revealed the strong association of rs443861 



with RA, indicating a genetic effect of TXNDC5 on RA 
risk. Although the genetic data of the present study 
indicated the possible association of TXNDC5 to RA, 
not enough data support the idea that the increased 
expression was caused by a genetic mechanism. The 
increased expression of TXNDC5 could be induced by 
hypoxia in RA rather than genetic variation of the gene. 
To determine whether variations in the TXNDC5 gene 
contributed to the risk of developing nonsegmental viti- 
ligo (NSV), Jeong et al, conducted a case-control asso- 
ciation study within a Korean population. They 
genotyped seven SNPs and found that three exonic 
SNPs (rsl043784, rs7764128 and rs8643) were statisti- 
cally associated with NSV. The haplotypes AGG and 
GAA, consisting of rsl043784, rs7764128 and rs8643, 
demonstrated a significant association with NSV [18]. 
Lin et al, reported that SNP rsl3873 and haplotypes 
rsl225934 to rsl3873 of BMP6-TXNDC5 genes were 
significantly associated with schizophrenia [19]. These 
reports indicate that TXNDC5 plays a role in the patho- 
genesis of other diseases. Our results demonstrated that 
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Figure 4 LD plot of the 96 genotyped SNPs in the TXNDC5 gene. (A) Linkage disequilibrium in tine RA group. (B) Linl<age disequilibrium in the AS 
group. Red areas representing higher levels of LD. Blue areas represent LD comparisons with low confidence of estimation. Dark triangles represent 
haplotype blocks. Numbers in squares are D' values. AS, ankylosing spondylitis; RA, rheumatoid arthritis; SNP, single nucleotide polymorphism. 



Table 4 Haplotype frequencies in a case-control cohort of patients with rheumatoid arthritis 



Haplotype 


Freq. 


Case, Control Ratio Counts 


Case, Control Frequencies 


P-Value 


Block 1 










AGG 


0.827 


426.9: 101.1, 273.0: 45.0 


0.808, 0.858 


0.0625 


GAA 


0.145 


81.0: 447.0, 42.0: 276.0 


0.153, 0.132 


0.3939 


AAG 


0.012 


7.8: 520.2, 2.0: 316.0 


0.015, 0.006 


0.2649 


Block 2 










AA 


0.814 


418.9: 109.1, 269.9: 48.1 


0.793, 0.849 


0.0446 


GG 


0.155 


87.9: 440.1, 42.9: 275.1 


0.166, 0.135 


0.2187 


AG 


0.024 


15.1: 512.9, 5.1: 312.9 


0.029, 0.016 


0.2441 


Block 3 










GAAG 


0.597 


297.7: 230.3, 207.0: 111.0 


0.564, 0.651 


0.0125 


AGGA 


0.236 


135.3: 392.7, 64.0: 254.0 


0.256, 0.201 


0.0674 


AGAG 


0.146 


83.8: 444.2, 40.0: 278.0 


0.159, 0.126 


0.1902 


AAAG 


0.012 


6.1: 521.9, 4.0: 314.0 


0.012, 0.013 


0.8831 
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Table 4 Haplotype frequencies in a case-control cohort of patients with rheumatoid arthritis (Continued) 



Block 4 
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block O 










A A 

AA 


0.71 5 
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n n "7ni 
U./23, U./U 1 


0.4774 


A r 


U. 1 o 


OO./. 43/.0, DO.J. zj4.j 
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U. 1 by, U.zUU 


U.zjb4 


GG 


0.1 


53.3: 472.7, 31.5: 286.5 


0.101, 0.099 


0.9173 


Block 7 










AA 


0.699 


374.8: 153.2, 216.9: 101.1 


0.710, 0.682 


0.3934 


AG 


0.214 


110.9: 417.1, 70.1: 247.9 


0.210, 0.220 


0.7211 


GG 


0.085 


41.1: 486.9, 30.9: 287.1 


0.078, 0.097 


0.3288 


Block 8 










AG 


0.531 


294.0: 232.0, 153.8: 164.2 


0.559, 0.484 


0.0336 


AA 


0.376 


187.8: 338.2, 129.2: 188.8 


0.357, 0.406 


0.1518 


GA 


0.092 


42.9: 483.1, 34.8: 283.2 


0.082, 0.109 


0.1742 


Table 5 Haplotype frequencies in a case-control cohort of patients with ankylosing spondylitis 


Haplotype 


Freq. 


Case, control ratio counts 


Case, control frequencies 


P-value 


Block 1 










AGGG 


0.755 


75.0: 27.0, 242.0: 76.0 


0.735, 0.761 


0.5994 


GAAA 


0.136 


15.0: 87.0, 42.0: 276.0 


0.147, 0.132 


0.7006 


AAGG 


0.1 


11.0: 91.0, 31.0: 287.0 


0.108, 0.097 


0.7616 


Block 2 










AAAGAAG 


0.541 


45.0: 57.0, 182.2: 135.8 


0.441, 0.573 


0.0198 


GAAAGGA 


0.23 


34.0: 68.0, 62.5: 255.5 


0.333, 0.196 


0.0043 


GGGAGAG 


0.126 


13.0: 89.0, 40.0: 278.0 


0.127, 0.126 


0.9686 


GAAGAAG 


0.065 


6.0: 96.0, 21.3: 296.7 


0.059, 0.067 


0.7731 


Block 3 










AGGAGGGGA 


0.609 


49.7: 51.9, 206.0: 112.0 


0.489, 0.648 


0.0044 


CGAGCAAAC 


0.187 


29.8: 71.8, 49.0: 269.0 


0.293, 0.154 


0.0018 


CAGGCAAGA 


0.134 


14.3: 87.3, 42.0: 276.0 


0.140, 0.132 


0.829 


CGAGCAAAA 


0.022 


1.2: 100.4, 8.0: 310.0 


0.012, 0.025 


0.4324 


CGAGCAAGA 


0.017 


4.0: 97.6, 3.0: 315.0 


0.039, 0.010 


0.0408 


CAAGCAAGA 


0.012 


0.0: 101.6, 5.0: 313.0 


0.000, 0.016 


0.2042 


Block 4 










GG 


0.664 


58.0: 44.0, 221.0: 97.0 


0.569, 0.695 


0.0187 


AA 


0.228 


28.0: 74.0, 67.7: 250.3 


0.275, 0.213 


0.1962 


AG 


0.108 


16.0: 86.0, 29.3: 288.7 


0.157, 0.092 


0.067 


Block 5 










AAA 


0.523 


45.9: 56.1, 173.7: 144.3 


0.450, 0.546 


0.0915 


AGC 


0.29 


33.3: 68.7, 88.6: 229.4 


0.327, 0.279 


0.3548 


GGC 


0.097 


7.6: 94.4, 33.0: 285.0 


0.075, 0.104 


0.3851 


AAC 


0.072 


10.7: 91.3, 19.4: 298.6 


0.105, 0.061 


0.1339 


AGA 


0.013 


4.1: 97.9, 1.3: 316.7 


0.040, 0.004 


0.0053 
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Table 5 Haplotype frequencies in a case-control cohort of patients with ankylosing spondylitis (Continued) 



Block 6 
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Block 8 
















0.49 


53.9: 


4o. 1 , 1 J 1 .o. 1 DD.Z 


0.529, 


0.477 


0.3656 


r' A A A A 

bAAAA 


0.1 98 
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oU.y, DZ. 1 . zbb.y 


0.206, 


0.1 95 
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A A A A 

AAAAb 
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0.1 4 


A rr' A A 
ALbAA 


0.057 


5.0: 


n"7 r> inn. inn n 

y/.u, ly.u. zyy.u 


0.049, 


0.060 


0.6831 


AAGAA 


0.042 


4.0: 


98.0, 13.9: 304.1 


0.039, 


0.044 


0.8369 


Block 9 














AT 


0.517 


53.0: 


49.0, 164.0: 154.0 


0.520, 


0.516 


0.9455 


GA 


0.28 


26.6: 


75.4, 91.0: 227.0 


0.261, 


0.286 


0.6165 


GT 


0.203 


22.4: 


79.6, 63.0: 255.0 


0.220, 


0.198 


0.6357 


Block 10 














AG 


0.802 


78.2: 


23.8, 258.5: 59.5 


0.767, 


0.813 


0.3066 


GA 


0.161 


18.7: 


83.3, 48.8: 269.2 


0.184, 


0.153 


0.4683 


GG 


0.032 


3.8: 


98.2, 9.5: 308.5 


0.037, 


0.030 


0.7032 



Table 6 Allele and genotype frequencies in a case-control cohort of patients with RA 


dbSNP identity 


Allele/ 
Genotype 


No. of patients with RA (%) 


No. of controls (%) 


Fisher's P-value 


rs2277105 


C 


1581 (0.840) 


1512 (0.848) 


0.491098 




O(freq) 


19 (0.010) 


12 (0.007) 






T 


283 (0.150) 


260 (0.146) 






CC 


670 (0.719) 


650 (0.734) 


0.549842 




CT 


241 (0.259) 


212 (0.239) 






JT 


21 (0.023) 


24 (0.027) 




rsl 1962800 


A 


1568 (0.831) 


1523 (0.853) 


0.173781 




O(freq) 


15 (0.008) 


11 (0.006) 






G 


304 (0.161) 


251 (0.141) 






AA 


656 (0.701) 


656 (0.740) 


0.184145 




AG 


256 (0.274) 


211 (0.238) 






GG 


24 (0.026) 


20 (0.023) 




rs443861 


A 


357 (0.190) 


275 (0.154) 


0.00832 




O(freq) 


19 (0.010) 


12 (0.007) 






G 


1507 (0.800) 


1497 (0.839) 






AA 


51 (0.055) 


26 (0.029) 


0.01011 




AG 


255 (0.274) 


223 (0.252) 






GG 


626 (0.672) 


637 (0.719) 




rs369086 


A 


485 (0.257) 


460 (0.258) 


0.510123 




O(freq) 


18 (0.010) 


11 (0.006) 






G 


1381 (0.733) 


1314 (0.736) 






AA 


61 (0.065) 


71 (0.080) 


0.260584 




AG 


363 (0.389) 


318 (0.359) 






GG 


509 (0.546) 


498 (0.561) 





898 controls and 951 RA cases observed. RA, rheumatoid arthritis; SNP, single nucleotide polymorphism 
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rsl043784, rs7764128, rsl225934 and rs8643 were not 
significantly associated with RA and AS. 

Conclusions 

Our study demonstrated significantly increased TXNDC5 
expression in the synovium and blood of RA patients, 
which may contribute to the irregular angiogenesis and 
abnormal cell differentiation observed in the synovial 
membrane. The study also revealed the genetic effect of 
TXNDC5 on RA and AS risk. 

Additional material 



Additional file 1: Supplementary materials and methods. This table 
summarizes the clinical data of patients with RA, OA and AS. 

Additional file 2: Supplementary results. This table provides the raw 
microarray data to perform association, LD and haplotype analysis. We 
genotyped 96 SNPs across the TXNDC5 gene from 267 Han Chinese 
patients with RA, 51 patients and 160 control individuals. All SNPs 
yielded genotype data, and the study sample success rate was 99.1%. 
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